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(57) ABSTRACT 

A method of determining kinetic parameters for a reversible 
molecular interaction between a ligand immobilized to a 

solid support surface and a binding partner to the ligand in 
solution, comprises sequentially, without intermediate 
regeneration or renewal of the immobilized ligand, flowing 
a plurality of fluid volumes containing different known 
concentrations of the binding partner over the solid support 
surface, monitoring the momentary amount of binding part- 
ner bound to the solid support surface related to time and 
solution concentration of binding partner and collecting the 
binding data, and determining the kinetic parameters by 
globally fitting a predetermined kinetic model for the inter- 
action between the binding partner and the immobilized 
ligand to the collected binding data, which model allows for 
mass transport limitation at the solid support surface. An 
analytical system for carrying out the method, a computer 
program, a computer program product and a computer 
system for performing the method are also disclosed. 
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METHOD AND SYSTEM FOR DETERMINATION 
OF MOLECULAR INTERACTION PARAMETERS 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

[0001] This application claims the benefit of U.S. Provi- 
sional Patent Application No. 60/477,909, filed Jun. 12, 
2003 and U.S. Provisional Patent Application No. 60/526, 
364, filed Dec. 1, 2003; and also claims priority to Swedish 
Patent Application No. 0301639-1, filed Jun. 6, 2003 and 
Swedish Patent Application No. 0303214-1, filed Dec. 1, 
2003; all of these applications are incorporated herein by 
reference in their entireties. 

BACKGROUND OF THE INVENTION 
[0002] 1. Field of the Invention 

[0003] The present invention relates to the determination 
of kinetic parameters for molecular interactions, and more 

particularly to a method for determining kinetic parameters 
for the interaction between a molecule immobilized to a 
solid support surface and a binding partner to the molecule 
in solution. The invention also relates to an analytical 
system, a computer program product and a computer system 
for performing the method. 

[0004] 2. Description of the Related Art 

[0005] Analytical sensor systems that can monitor inter- 
actions between molecules, such as biomolecules, in real 
time are gaining increasing interest. These sensor systems, 
usually referred to as interaction analysis sensor systems or 
biospecific interaction analysis sensor systems, are often 
based on optical biosensors and affinity analysis and offer a 
rapid way to determine in real time inter alia equilibrium and 
rate constants without the need to label the interacting 
molecules. They have been used in the study of a variety of 
biomolecules, including proteins, nucleic acids, lipids and 
carbohydrates. In these systems, a sensor surface having one 
of the molecular reactants immobilized thereto is contacted 
with a solution containing the other reactant, either by 
providing a flow of the solution past the sensor surface, or 
in a cuvette or the like, and binding interactions at the 
surface are detected. 

[0006] Conventionally, to determine, for example, asso- 
ciation and dissociation rate constants (k^ and k^, respec- 
tively) for the interaction between two interacting mol- 
ecules, one of the molecules, often referred to as the ligand, 
is immobilized to a sensor surface and the other molecule, 
often referred to as the analyte, is provided in solution at 
several different known concentrations. Each concentration, 
or sample, of the analyte is then contacted with the sensor 
surface, either in a laminar flow past the sensor surface, or 
in a cuvette or the like, to permit association of the analyte 
to the sensor surface. After a sample has been brought to 
contact the sensor surface, the surface is contacted with a 
solution free from analyte, usually buffer, to permit disso- 
ciation of the analyte from the immobilized ligand. During 
these association and dissociation phases, the amount of 
binding of analyte to the surface is continuously detected 
and the binding data is collected. Before contacting the 
sensor surface with sample of a new analyte concentration, 
the ligand surface is restored or "regenerated" by treating the 
surface with a regeneration solution capable of removing 



any bound analyte while not destroying the ligand. In that 
way, all the different samples will contact essentially one 

and the same ligand surface as far as ligand density is 
concerned. The association and dissociation rate constants 
can then be obtained from the collected binding data by 
fitting the data to mathematical descriptions of interaction 
models in the form of differential equations. Usually, the 
binding data for all the samples are used in the same fit, a 
procedure referred to as global fitting. From the determined 
association and dissociation rate constants k^ and k^, the 
equilibrium constant, Kj^, and the affinity constant 
(Ka=1/Kj5) of the interaction can in turn be calculated. 
Alternatively, provided that the interaction reaches steady 
state during the association phase, the equilibrium constant 
can be directly obtained from the binding data without 
fitting. 

[0007] Problems may arise, however, when the ligands are 
covalently immobilized to the sensor surface and suitable 
regeneration conditions are difficult to find. Renewed bind- 
ing of the ligand via an immobilized capture agent before the 
contacting with each new analyte concentration could then 
be an alternative, but has the disadvantage of consuming 
large amounts of ligand for the determination. 

[0008] Determination of equilibrium constants by a titra- 
tion procedure without the requirement for regeneration of 
the immobilized ligands is described by Schuck, P., et al. 
(199S) Anal Biochem. 265, 79-91. The sample is continu- 
ously circulated in a closed loop over two sensor spots of a 
commercial surface plasmon resonance biosensor. One of 
the sensor spots is functionalized with an immobilized 
ligand for a soluble analyte in the sample, and the other 
sensor spot serves as a reference surface. A binding isotherm 
for the interacting molecules is obtained by stepwise equi- 
librium titration of the analyte into the circulating loop, i.e., 
the sensor spots are sequentially contacted with stepwise 
increased concentrations of analyte until equilibrium is 
attained for each concentration. This equilibrium titration is 
said to be especially useful for the determination of binding 
constants in high-affinity systems since it eliminates the 
need for interpretation of binding kinetics and thereby 
problems that may arise from mass transport limitations. 

[0009] A similar stepwise equilibrium titration procedure 

is described for a cuvette -based biosensor design in Hall, D. 
R., and Winzor, D. J. {1991) Anal Biochem. 244, 152-160. 

[0010] Also Myszka, D. G., et al. (1998) Anal. Biochem. 
265, 326-330 discloses equilibrium analysis of high affinity 
interactions using a surface plasmon resonance-based bio- 
sensor. In this approach, the time available to collect asso- 
ciation phase data is increased by placing the analyte 
directly into the running buffer. Complete equilibrium bind- 
ing profiles were generated without a regeneration step by 
changing the concentration of analyte and allowing the 
surface reactions to reequilibrate. Analyte concentrations 
were also decreased to demonstrate that the binding reac- 
tions were fully reversible. In this way, equilibrium disso- 
ciation constants for very high affinity interactions could be 
determined. 

[0011] Shank-Retzlaff, M. L., and SHgar, S. G. (2000) 
Anal Chem. 72, 4212-4220 describes a one-step method for 
determining kinetic rates and equilibrium binding affinities 
by a technique called analyte-gradient surface plasmon 
resonance SPR (AG-SPR) which eliminates the need for 
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surface regeneration. A gradient of analyte is passed over the 
sensor surface under continuous-flow conditions so that the 

concentration of analyte increases linearly with time. The 
rate at which analyte binds to immobilized ligands on the 
sensor surface is monitored by monitoring the change in the 
surface plasmon resonance as the analyte concentration 
increases. Kinetic rates are determined by fitting the data to 
a two-compartment model for the molecular interaction 
which permits use also for systems under mass transport 
limitations. 

[0012] Titration procedures to determine binding capacity 
and regeneration conditions have been proposed by Karls- 
son, R., et al., Biasymposium 1998, Edinburgh 2-4 Septem- 
ber 1998, 5th Biasymposium in Japan, Tokyo, Nov. 5-6, 
1998. 

[0013] US-Al-2003/014365 discloses a method for deter- 
mining interaction parameters, including rate constants, 
between an analyte and a Hgand immobilized to a sensor 
surface, where a measurement can be performed several 
times in succession, e.g., in a cuvette, with a stepwise 
modification of the analyte concentration each time. The 
measurements need not be completed to equilibrium but can 
be interrupted earlier and the concentration of the analyte 
can be raised or lowered. Separate fitting of the binding 
curve part corresponding to each analyte concentration is 
used to determine respective initial binding rates, from 
which the association and dissociation rate constants for the 
interaction are then determined. 

[0014] Thus, while titration procedures in combination 
with biosensors and affinity analysis are known per se in the 
art for determining equilibrium binding affinities, the use of 
multiple titrations to determine kinetic rate constants seems 
to be disclosed only in the above-mentioned publications 
Shank-Retzlaff, M. L., and SHgar, S. G. (2000), and US-Al- 
2003/014365. The method according to Shank-Retzlaff, M. 
L., and Sligar, S. G. (2000) suffers from the limitation that 
it requires the use of a continuous gradient of the analyte. 
The method of US-Al-2003/014365, on the other hand, is 
disadvantageous in that initial binding rates are frequently 
lower than the kinetic binding rate due to transport limita- 
tions, making initial binding rates unreliable for kinetic 
analysis. Further, this type of evaluation restricts the injec- 
tion order of analyte concentrations. 

[0015] From the prior art it may therefore be concluded 
that for determining kinetic rates for molecular interactions 
using systems based on biosensors and affinity analysis, it is 
necessary to regenerate the immobilized ligand prior to 
contacting the sensor surface with a different concentration 
of analyte to thereby present essentially one and the same 
ligand surface to each analyte concentration, unless (i) a 
continuous gradient of the analyte is used, or (ii) initial 
binding rates are determined in systems free from mass 
transport limitations. 

[0016] It is an object of the present invention to provide a 
sensor-based method for determining chemical interaction 
parameters, including kinetic rate constants, by stepwise 
titration, which method obviates regeneration procedures 
while permitting measurements under mass transport limi- 
tation. 

BRIEF SUMMARY OF THE INVENTION 

[0017] According to the present invention, it has surpris- 
ingly been found that the above and other objects and 



advantages can be achieved by an affinity sensor-based 
method for determining chemical interaction parameters, 
which method comprises performing, in one and the same 
experimental cycle, stepwise changes of the analyte concen- 
tration without intermediate regeneration of the sensor sur- 
face, and determining the interaction parameters using glo- 
bal analysis of the whole set of detected binding data with 
regard to a kinetic model for the molecular interaction that 
allows for mass transport limitation, such as, e.g., the 
afore-mentioned two-compartment model. In this method, 
kinetic data can thus be obtained in spite of mass transport 
hmitations and the different analyte concentrations may be 
used in any order. 

[0018] Therefore, in one aspect, the present invention 
provides a method of determining kinetic parameters for a 
reversible molecular interaction between a ligand immobi- 
lized to a sensor surface and a binding partner (analj^e) to 
the ligand in solution, comprising the steps of: 

[0019] a) sequentially, without intermediate regen- 
eration or renewal of the immobilized ligand, fiow- 
ing a plurality of fluid volumes containing different 
known concentrations of the binding partner over the 
solid support surface to permit association of binding 
partner to the immobilized ligand, 

[0020] b) flowing over the solid support surface a 
fluid volume free from binding partner to permit 
dissociation of binding partner from the ligand, 

[0021] c) monitoring during steps a) and b) the 
momentary amount of binding partner bound to the 
solid support surface related to time and solution 
concentration of binding partner and collecting the 
binding data, and 

[0022] d) determining the kinetic parameters by fit- 
ting, preferably globally, a predetermined kinetic 
model for the interaction between the binding partner 

and the immobilized ligand to the collected binding 
data, which model allows for mass transport limita- 
tion at the solid support surface. 

[0023] In a variant method, the different concentrations of 
binding partner are obtained by generating a gradient of the 
binding partner and intersecting segments of the gradient 
with segments of a fluid free from binding partner. 

[0024] In another aspect, the present invention therefore 
provides a method of determining kinetic parameters for a 
reversible molecular interaction between a ligand immobi- 
lized to a sensor surface and a binding partner (analyte) to 
the ligand in solution, comprising the steps of: 

[0025] a) sequentially, without intermediate regen- 
eration or renewal of the immobilized ligand, flow- 
ing a plurality of fluid volumes containing different 
concentrations of the binding partner over the solid 
support surface to permit association of binding 
partner to the immobilized ligand, wherein the fluid 
volumes containing binding partner are discrete seg- 
ments of a concentration gradient of the binding 
partner separated by segments of fluid free from 
binding partner, 

[0026] b) monitoring during step a) the momentary 
amount of binding partner boimd to the solid support 
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surface related to time and solution concentration of 
binding partner and collecting the binding data, and 

[0027] c) determining the kinetic parameters by fit- 
ting a predetermined kinetic model for the interac- 
tion between the binding partner and the immobi- 
lized ligand to the collected binding data. 

[0028] Preferably, the model allows for mass transport 
limitation at the solid support surface. 

[0029] In the case that the concentrations of all the gra- 
dient segments are not known, the unknown concentrations 
may be estimated. This may, for example, be done by having 
the fitting in step d) include local fitting of concentrations. 

[0030] The solid support surface is preferably a sensor 
surface, i.e., one which produces a detectable signal in 
response to a binding interaction at the surface. 

[0031] The term "monitoring" as used herein means that 
detection is performed at at least a plurality of times during 
steps a) and b), preferably at a large number of times. 

[0032] In another aspect, the present invention provides a 
method for determining the concentrations of binding part- 
ner in a plurality of binding partner-containing samples. 

[0033] In still another aspect, the present invention pro- 
vides a method for determining kinetic parameters for a 
plurality of different binding partners. 

[0034] In another aspect, the present invention provides an 
analytical system for studying molecular interactions, which 
system comprises data processing means for performing at 
least one of the methods. 

[0035] In yet another aspect, the present invention pro- 
vides a computer program comprising program code means 
for performing at least one of the methods. 

[0036] In still another aspect, the present invention pro- 
vides a computer program product comprising program code 
means stored on a computer readable medium or carried on 
an electrical or optical signal for performing at least one of 
the methods. 

[0037] In another aspect, the present invention provides a 
computer system containing a computer program compris- 
ing program code means for performing at least one of the 

methods. 

[0038] Other advantages, novel features and objects of the 
invention will become apparent from the following detailed 
description of the invention when considered in conjunction 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0039] FIG. 1 is a schematic side view of a biosensor 

system based on surface plasmon resonance (SPR). 

[0040] FIG. 2 is a representative sensorgram with a bind- 
ing curve having association and dissociation phases. 

[0041] FIG. 3 is a representative sensorgram obtained by 
sequential injections of a kinase inhibitor (staurosporin) over 
a surface with immobilized kinase according to the method 

of the present invention. 

[0042] FIG. 4 is the sensorgram of FIG. 3 having overlaid 
thereon the binding curve resulting from global fitting. 



[0043] FIG. 5 is a schematic illustration of a system for 
creating a pulsed concentration gradient. 

[0044] FIG. 6 is a schematic illustration of an alternative 

system for creating a pulsed concentration gradient. 

[0045] FIG. 7 is an overlay plot of observed and fitted 
binding data for sequential injections of a carbonic anhy- 
drase inhibitor (acetazolemide) over a surface with immo- 
bilized carbonic anhydrase. 

[0046] FIG. 8 is a similar overlay plot to that of FIG. 7 for 
sequential injections of another carbonic anhydrase inhibitor 
(azosulfamide). 

[0047] FIG. 9 is a similar overlay plot to those of FIGS. 
7 and 8 for sequential injections of still another carbonic 
anhydrase inhibitor (benzenesulfonamide). 

[0048] FIG. 10 is an overlay plot of observed and fitted 
binding data for sequential injections of a thrombin inhibitor 
(melagatran) over a surface with immobilized thrombin. 

[0049] FIG. 11 is an overlay plot of four binding curves 
and fitted data, three curves (to the left) representing sen- 
sorgrams with single injections of analyte, and one curve (to 
the right) representing a sensorgram for sequential injections 
of analyte. 

[0050] FIG. 12 is a sensorgram obtained by pulse injec- 
tion of camel antibody (SGS) over a surface with immobi- 
lized lysozyme. 

DETAILED DESCRIPTION OF THE 
INVENTION 

[0051] Unless defined otherwise, all technical and scien- 
tific terms used herein have the same meaning as commonly 
understood by a person skilled in the art related to this 
invention. Also, the singular forms "a", "an", and "the" are 
meant to include plural reference unless it is stated other- 
wise. 

[0052] As mentioned above, the present invention relates 
to the determination of molecular interaction parameters, 
including kinetic rate constants, for the interaction between 
a molecule immobilized to a solid support surface and a 
binding partner to the molecule in solution by a novel 
titration type method, preferably in combination with sensor 
based technology to study the molecular interactions and 
present the results in real time, as the interactions progress. 
Before describing the present invention in more detail, 
however, the general context in which the invention is used 
will be described. 

[0053] Chemical sensors or biosensors are typically based 
on label-free techniques, detecting a change in a property of 
a sensor surface, such as, e.g., mass, refractive index, or 
thickness for the immobilized layer, but there are also 
sensors relying on some kind of labelling. Typical sensor 
detection techniques include, but are not limited to, mass 
detection methods, such as optical, ther mo -optical and 
piezoelectric or acoustic wave, (including, e.g., surface 
acoustic wave (SAW) and quartz crystal microbalance 
(QCM)) methods, and electrochemical methods, such as 
potentiometric, conductometric, amperometric and capaci- 
tance/impedance methods. With regard to optical detection 
methods, representative methods include those that detect 
mass surface concentration, such as reflection-optical meth- 
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ods, including both external and internal reflection methods, 
angle, wavelength, polarization, or phase resolved, for 
example evanescent wave ellipsometry and evanescent 
wave spectroscopy (EWS, or Internal Reflection Spectros- 
copy), both may include evanescent field enhancement via 
surface plasmon resonance (SPR), Brewster angle refracto- 
metry, critical angle refractometry, frustrated total reflection 
(FTR), scattered total internal reflection (STIR) — ^which 
may include scatter enhancing labels, optical wave guide 
sensors; external reflection imaging, evanescent wave-based 
imaging such as critical angle resolved imaging, Brewster 
angle resolved imaging, SPR-angle resolved imaging, and 
the like. Further, photometric and imaging/microscopy 
methods, "per se" or combined with reflection methods, 
based on for example surface enhanced Raman spectroscopy 
(SERS), surface enhanced resonance Raman spectroscopy 
(SERRS), evanescent wave fluorescence (TIRF) and phos- 
phorescence may be mentioned, as well as waveguide inter- 
ferometers, waveguide leaking mode spectroscopy, reflec- 
tive interference spectroscopy (RIfS), transmission 
interferometry, holographic spectroscopy, and atomic force 
microscopy (AFR). 

[0054] Commercially available biosensors include the 
BIACORE® system instruments, manufactured and mar- 
keted by Biacore AB, Uppsala, Sweden, which are based on 
surface plasmon resonance (SPR) and permit monitoring of 
surface binding interactions in real time between a bound 
ligand and an analyte of interest. 

[0055] While in the detailed description and Examples that 
follow, the present invention is illustrated in the context of 
SPR spectroscopy, and more particularly the BIACORE® 
system, it is to be understood that the present invention is not 
limited to this detection method. Rather, any affinity-based 
detection method where an analyte binds to a ligand immo- 
bilized on a sensing surface may be employed, provided that 
a change at the sensing surface can be measured which is 
quantitatively indicative of binding of the analyte to the 
immobilized ligand thereon. 

[0056] The phenomenon of SPR is well known, suffice it 
to say that SPR arises when light is reflected under certain 
conditions at the interface between two media of different 
refractive indices, and the interface is coated by a metal film, 
typically silver or gold. In the BIACORE® instruments, the 
media are the sample and the glass of a sensor chip which 
is contacted with the sample by a microfluidic flow system. 
The metal film is a thin layer of gold on the chip surface. 
SPR causes a reduction in the intensity of the reflected light 
at a specific angle of reflection. This angle of minimum 
reflected light intensity varies with the refractive index close 
to the surface on the side opposite from the reflected light, 
in the BIACORE® system the sample side. 

[0057] A schematic illustration of the BIACORE® system 
is shown in FIG. 1. Sensor chip 1 has a gold film 2 
supporting capturing molecules 3, e.g., antibodies, exposed 
to a sample flow with analytes 4 (e.g., an antigen) through 
a flow channel 5. Monochromatic p-polarised light 6 from a 
light source 7 (LED) is coupled by a prism 8 to the 
glass/metal interface 9 where the light is totally reflected. 
The intensity of the reflected light beam 10 is detected by an 
optical detection unit (photodetector array) 11. 

[0058] A detailed discussion of the technical aspects of the 
BIACORE® instruments and the phenomenon of SPR may 



be found in U.S. Pat. No. 5,313,264. More detailed infor- 
mation on matrix coatings for biosensor sensing surfaces is 

given in, for example, U.S. Pat. Nos. 5,242,828 and 5,436, 
161. In addition, a detailed discussion of the technical 
aspects of the biosensor chips used in connection with the 
BIACORE® instruments may be found in U.S. Pat. No. 
5,492,840. The fuU disclosures of the above-mentioned U.S. 
patents are incorporated by reference herein. 

[0059] When molecules in the sample bind to the captur- 
ing molecules on the sensor chip surface, the concentration, 
and therefore the refractive index at the surface changes and 
an SPR response is detected. Plotting the response against 
time during the course of an interaction will provide a 
quantitative measure of the progress of the interaction. Such 
a plot is usuafly called a sensorgram. In the BIACORE® 
system, the SPR response values are expressed in resonance 
units (RU). One RU represents a change of 0.0001° in the 
angle of minimum reflected light intensity, which for most 
proteins and other biomolecules correspond to a change in 
concentration of about 1 pg/mm^ on the sensor surface. As 
sample containing an analj^e contacts the sensor surface, the 
ligand bound to the sensor surface interacts with the analyte 
in a step referred to as "association." This step is indicated 
on the sensorgram by an increase in RU as the sample is 
initially brought into contact with the sensor surface. Con- 
versely, "dissociation" normally occurs when the sample 
flow is replaced by, for example, a bufter flow. This step is 
indicated on the sensorgram by a drop in RU over time as 
analyte dissociates from the surface-bound ligand. 

[0060] A representative sensorgram (binding curve) for a 
reversible interaction at the sensor chip surface is presented 
in FIG. 2, the sensing surface having an immobilized 
capturing molecule, or ligand, for example an antibody, 
interacting with a binding partner therefor, or analyte, in a 
sample. (The detection curves, or sensorgrams, produced by 
biosensor systems based on other detection principles men- 
tioned above will have a similar appearance.) The y-axis 
indicates the response (here in resonance units, RU) and the 
X-axis indicates the time (here in seconds). Initially, buffer is 
passed over the sensing surface giving the baseline response 
A in the sensorgram. During sample injection, an increase in 
signal is observed due to binding of the analyte. This part B 
of the binding curve is usually referred to as the "association 
phase". EventuaUy, a steady state condition is reached where 
the resonance signal plateaus at C (this state may, however, 
not always be achieved). At the end of sample injection, the 
sample is replaced with a continuous flow of buffer and a 
decrease in signal reflects the dissociation, or release, of 
analyte from the surface. This part D of the binding curve is 
usually referred to as the "dissociation phase". The analysis 
is ended by a regeneration step where a solution capable of 
removing bound analyte from the surface, while (ideally) 
maintaining the activity of the ligand, is injected over the 
sensor surface. This is indicated in part E of the sensorgram. 
Injection of buffer restores the baseline A and the surface is 
now ready for a new analysis. 

[0061] From the profiles of the association and dissocia- 
tion phases B and D, repectively, information regarding the 
binding and dissociation kinetics is obtained, and the height 
of the resonance signal represents surface concentration 
(i.e., the response resulting from an interaction is related to 
the change in mass concentration on the surface). This will 
now be explained in more detafl below. 
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[0062] Assume a reversible reaction between an analyte A 
and a surface -bound (immobilized) capturing molecule, or 
ligand, B which is not diffusion or mass transfer limited and 
obeys pseudo first order kinetics: 

A+BoAB 

[0063] This interaction model (usually referred to as the 
Langmuir model), which assumes that the analyte (A) is 
both monovalent and homogenous, that the ligand (B) is 
homogenous, and that all binding events are independent, is 
in fact applicable in the vast majority of cases. 

[0064] The rate of change in surface concentration of 
analyte A (=rate of change in concentration of formed 
complex AB) during analyte injection is the sum of the rates 
of the analyte A going on and off: 



k^C+kd and the vertical intercept is k^RmaxC. If the bulk 
concentration C is known and R^^ax has been determined 
(e.g., by saturating the surface with a large excess of 
analyte), the association rate constant k^^ and the dissociation 
rate constant k^ can be calculated. A more convenient 
method is, however, fitting of the integrated function (5), or 
numerical calculation and fitting of the differential Equation 
(4), preferably by means of a computer program as will be 
described below. 

[0071] The rate of dissociation can be expressed as: 



d[AB] 

~~dr '' 



(1) 



[0065] where [A] is the concentration of analyte A, [B] is 
the concentration of the ligand B, [AB] is the concentration 
of the reaction complex AB, k^ is the association rate 
constant, and k^ is the dissociation rate constant. 

[0066] After a time t, the concentration of unbound ligand 
B at the surface is [BT]-[AB], where [BT] is the total, or 
maximum, concentration of ligand B. Insertion into Equa- 
tion (1) gives: 



d[AB] 

~dr 



kAA]{[BT]-[AB]}-kAAB] 



(2) 



[0067] In terms of detector response units (AB is 
detected), this can be expressed as 



dR 

= kaC{Rjnax. - R) - kdR 



(3) 



[0068] where R is the response at time t in resonance units 
(RU), C is the initial, or bulk, concentration of free analyte 
(A) in solution, and R^ax the response (in RU) obtained 
if analyte (A) had bound to all ligand (B) on the surface. 
Rearrangement of Equation (3) gives: 



[0072] and in integrated form: 



(7) 



[0073] where Rq is the response at the beginning of the 
dissociation phase (when the buffer wash of the surface 
starts). 

[0074] Equation (6) may be linearized: 



-kdt 



(8) 



[0075] and a plot of ln[R/Ro] versus t will produce a 
straight line with the slope^-k^. More conveniently, how- 
ever, the dissociation rate constant k^ is determined by fitting 

the exponential rate equation (7). 

[0076] Affinity is expressed by the association constant 
K^^k^/kd, or the dissociation constant (also referred to as the 

equilibrium constant) Kj5=kd/ka. 

[0077] The association constant may alternatively be 
determined from Equation (3), where dR/dt=0 at equilib- 
rium, giving: 

A:Aq=^aC(i?n.ax-^eq) (9) 

[0078] where R^^ is the detector response at equilibrium. 
Since k^/kd^K^, insertion in Equation (9) and rearrangement 

gives: 



— = -K^Rec^ 



(10) 



=kaCR^-{kaC-\-kd)R 



[0069] where R is the response in resonance units (RU). In 
integrated form, the equation is: 



[0070] Now, according to equation (4), if dR/dt is plotted 
against the bound analyte concentration R, the slope is 



[0079] If binding reactions are performed at multiple 

concentrations, the data may either be fitted or Rgq/C may be 
plotted versus R^^ which gives the slope =-K^. Such an 
equilibrium analysis may be performed when the rates of 
association and dissociation are too rapid to measure accu- 
rately. 

[0080] To obtain reliable kinetic constants, the above 
described analysis is usually repeated for a number of 
different analyte concentrations and, suitably, also at least 
one other ligand density at the sensor surface. 

[0081] Software for the analysis of kinetic and other 
biosensor data is commercially available. Thus, for example. 
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evaluation of kinetic data produced by the BIACORE® 
instruments is usually performed with the dedicated BIAe- 
valuation software (supplied by Biacore AB, Uppsala, Swe- 
den) using numerical integration to calculate the differential 
rate equations and non-linear regression to fit the kinetic 
parameters by finding values for the variables that give the 
closest fit reducing the sum of squared residuals to a 
minimum. The residuals are the difference between the 
calculated and the experimental curve at each point, squared 
residuals being used to weight equally deviations above and 
below the experimental curve. The simi of squared residual 
is expressed by Equation (11): 



(11) 



[0082] where S is the sum of squared residuals, r^ is the 
fitted value at a given point, and r^ is the experimental value 
at the same point. 

[0083] For example, for the molecular interaction 
described above, such software -assisted data analysis is 
performed by, after subtracting background noises, making 
an attempt to fit the above-mentioned simple 1:1 Langmuir 
binding model as expressed by Equations (5) and (7) above 
to the measurement data. 

[0084] Usually the binding model is fitted simultaneously 
to multiple binding curves obtained with different analyte 

concentrations C (and/or with different levels of surface 
derivatization R^ax)- This is referred to as "global fitting", 
and based on the sensorgram data such global fitting estab- 
lishes whether a single global k^ or k^ will provide a good fit 
to all the data. The results of the completed fit are presented 
to the operator graphically, displaying the fitted curves 
overlaid on the original sensorgram curves. The closeness of 
the fit is also presented by the chi-squared (x^) value, a 
standard statistical measure. For a good fitting, the chi- 
squared value is in the same magnitude as the noise in RU^. 
Optionally, "residual plots" are also provided which give a 
graphical indication of how the experimental data deviate 
from the fitted curve showing the difference between the 
experimental and fitted data for each curve. The operator 
then decides if the fit is good enough. If not, the sensorgram 
or sensorgrams exhibiting the poorest fit are excluded and 
the fitting procedure is run again with the reduced set of 
sensorgrams. This procedure is repeated until the fit is 
satisfactory. 

[0085] Sometimes, the above-mentioned 1:1 binding reac- 
tion model will not be valid, which requires the data set to 

be reanalysed using one or more other reaction models. Such 
alternative models may include, for example, a one to one 
reaction influenced by mass transfer (mass transport), two 
parallel independent one to one reactions, two competing 
reactions, and a two state reaction. Parallel reactions can 
occur when the immobiHzed ligand is heterogeneous, 
whereas a heterogenous analyte may give rise to competing 
reactions. A two state reaction indicates a conformation 
change that gradually leads to a more stable complex 
between ligand and analyte. For differential rate equations 
reflecting these alternative reaction models, it may be 
referred to, for example, Karlsson, R., and Fait, A. (1997) J. 



Immunol Methods 200, 121-133 (the disclosure of which is 
incorporated by reference herein). For a more comprehen- 
sive description of curve fitting with regard to the BIA- 
CORE® system, it may be referred to the BIAevaluation 
Software Handbook (Biacore AB, Uppsala, Sweden) (the 
disclosure of which is incorporated by reference herein). 

[0086] As regards mass transport, transport effects will 
influence the kinetics of binding when the reaction rate is 
fast compared to the rate of transport. For a reaction where 
analyte (A) binds to immobilized ligand (B), binding with 
mass transport Hmitation may be represented by the reaction 
formula: 

[0087] The differential equations describing the binding 
interaction will therefore include terms for mass transfer of 
analyte to the surface. For a flow cell, a "two -compartment" 
model consisting of a set of coupled ordinary differential 
equations and described in inter alia Myszka, D. G., et al. 
(1998) Biophys. J. 75, 583-594, and Shank-Retzlaff, M. L., 
and Sligar, S. G. (2000) supra (the relevant disclosures of 
which are incorporated by reference herein) is considered to 
give a reasonable description of the binding kinetics when 
the data are influenced by mass transport. In this model, the 
flow cell is assumed to be divided into two compartments, 
one in which the concentration of analyte is constant, and a 
second near the sensor surface where the analyte concen- 
tration depends on the mass transport rate, the siu"face 
density of ligand, and the reaction rate constants. 

[0088] For the interaction of a monovalent analyte (A) 
reacting with an immobilized monovalent ligand (B), this 
model may be represented by the following two differential 
equations (12) and (13): 



dA 
It 

dAB 
"dT '' 



— = i-kaA(BT - AB) + kdAB + kt(C - A)) 
dt 



■ kaAiBr - AB) - k^AB 



(12) 
(13) 



[0089] where k^ is the transport coefficient describing 
diffusive movement of analyte between the compartments, 
B^ is the total ligand concentration, A is the bulk concentra- 
tion of free analyte, C is the injection (i.e., initial) analyte 
concentration, AB is the concentration of complex AB 
(=surface density of bound analyte), and k^ and k^ are the 
association and dissociation rate constants, respectively. 

[0090] It is understood that the use of this two-compart- 
ment model analyzing the kinetic globally will extend the 
range of reaction rates that can be determined with instru- 
ments like the BIACORE® and has also been incorporated 
in the above-mentioned BIAevaluation analysis software 
(available from Biacore AB, Uppsala, Sweden). 

[0091] Theoretical interaction data for different interaction 
models may be simulated using a simulator, such as, e.g., the 
simulator which is built-in in the above-mentioned BIAe- 
valuation analysis software. In this way, the effects of 
varying experimental parameters may be explored, thereby 
helping to determine which interaction conditions that may 
be crucial in distinguishing between alternative kinetic mod- 
els. 
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[0092] The invention 

[0093] As mentioned above, determining kinetic param- 
eters using, for example, a BIACORE® instrument, con- 
ventionally comprises monitoring the interaction of analyte 
with immobilized ligand for a plurality of different analyte 
concentrations with regeneration, or renewal, of the immo- 
bilized ligand before each measurement with a new analyte 
concentration. A kinetic model for the interaction is then 
fitted globally to the collected binding data, usually in the 
form of sensorgrams, to give the kinetic parameters. 

[0094] According to the present invention, such a flow 
cell-based analysis for determining kinetic parameters may 
be substantially improved and speeded up by a titration type, 
or "sequential injection", procedure wherein the ligand- 
supporting surface is successively contacted, in one and the 
same analytical cycle, with the different analyte concentra- 
tions to produce a continuous sensorgram. A kinetic inter- 
action model allowing for mass transport limititation, such 
as the two -compartment model mentioned above, is then 
fitted globally to the entire sensorgram to calculate the 
kinetic parameters. In addition to being considerably less 
time-consuming, this titration type method reduces the 
amount of experimental data to be evaluated and eliminates 
the risk of regeneration conditions destroying immobilized 
ligand. 

[0095] The order of introducing the different analyte con- 
centrations is not critical to the success of the inventive 
procedure. Rather, the analyte concentrations may, in fact, 
be introduced in any order. For example, the analyte con- 
centrations may be successively increased, or successively 
decreased, or may alternately be higher and lower, etc. 
Further, the same concentration of analyte may be intro- 
duced repeatedly, if desired. 

[0096] The injections of analyte, which preferably may be 
short, for a BIACORE® or similar instrument, e.g., in the 
order of 30-60 seconds, may be interrupted when an 
increased analyte concentration only leads to a marginally 
increased response. This may advantageously be effected 
using so-called adaptive software. The formed complex is 
then preferably allowed to dissociate for a longer period of 
time, e.g., from about 5 to 60 minutes in the context of a 
BIACORE® or similar instrument, to add adequate disso- 
ciation data for the kinetic evaluation. Such a dissociation 
period may, optionally, instead (or additionally) be per- 
formed between any two analyte injections in the cycle. 

[0097] It is to be noted that the design of the above- 
mentioned BIACORE® instruments, for example, requires 
that each injection of analyte be followed by injection of 
nmning buffer, permitting partial dissociation of the bound 
analyte, before the next concentration of analyte is injected. 
It is understood that such alternating analyte-free liquid 
volumes, for example with about one minute intervals, will 
add dissociation data to the those obtained by the above- 
mentioned (preferably longer) dissociation period. Provided 
that sufficient dissociation data is obtained by these alter- 
nating injections of analyte-free fluid, the longer dissocia- 
tion period may even be dispensed with. Intermediate bufler 
injections may, of course, also be used, even if a sensor 
instrument design as such permits sequential injections of 
analyte without buffer injections between them. 
[0098] At least for some applications, it is preferred that 
the flow is substantially constant during the entire experi- 
mental cycle. 



[0099] The sequential injection procedure may optionally 
be repeated at least once with the same or another sequence 

of analyte concentrations, and all the collected binding data 
may then be fitted simultaneously to the kinetic interaction 
model. 

[0100] Optionally, the fluid volumes containing the differ- 
ent analyte concentrations may be different volume seg- 
ments of a concentration gradient of the analyte. A method 
for the sequential injection of segments of a decreasing (or 
increasing) analyte gradient is described in our copending 
U.S. application entitled "Method and apparatus for char- 
acterization of interactions" (the disclosure of which is 
incorporated by reference herein). 

[0101] In this "pulse injection" method, each injection 
consists of a series of short sample pulses, suitably 4 or 5 up 
to 40 pulses, preferably 15-30, more preferably about 20 
pulses, generated by alternating flows of analyte-containing 
sample, and another liquid, such as buffer. Each pulse 
preferably has a volume of 1-40 //I, preferably 10-40 //I, 
more preferably 15-25 ju\, suitably about 20 fd. The duration 
of the pulses, i.e., each segment of solution can be 8-20, 
preferably about 10-15, suitably 12 seconds long, and the 
flow rate for the sample liquid through the flow cell may be 
50-200, preferably 80-120, suitably 100 /il/min. 

[0102] In this way, information from several concentration 
levels will be generated through a single injection in that 
each pulse of the injection in principle constitutes one 
concentration. Optionally, some pulses during an injection 
may be discarded, whereby the discarded segments will not 
be passed through the sensor. Alternatively, some aliquot(s) 
of liquid can be discarded even before performing the 
alternating buffer injections to create the separated seg- 
ments. 

[0103] Depending on how the concentration gradient is 
prepared, the concentrations of all segments may or may not 
be known. If all concentrations are known, the kinetic 
parameters may be determined by global fitting as described 
above. If, on the other hand, only the concentration of, say, 
an edge segment(s) is known, such as in the case of a 
"dispersion gradient", the fitting may have to include local 
fitting of the analyte concentrations. 

[0104] If desired, binding data from one or more sequen- 
tial experiments performed as described above may be 
combined with binding data from one or more conventional 
type single concentration experiments, and the combined 

data are then fitted to the kinetic interaction model. Thus, 
generally, the total data set that is globally fitted to the 
kinetic interaction model may comprise any combination of 
binding curves including one or more "sequential" binding 
curves and one or more single concentration binding curves. 

[0105] To facilitate the understanding of the invention, 
reference is now made to FIG. 3 which shows an exemplary 
sensorgram (as obtained with the above-mentioned BIA- 
CORE® instrument) for a titration with successively 
increased concentrations of an analyte (A) reacting with an 
immobflized ligand (B) to form a complex (AB) which gives 
rise to the detected SPR response. This titration is described 
in more detail in Example 1 below. Referring to the sensor- 
gram, analyte concentration 80 nM was injected between 
404-464 seconds, 240 nM between 523-583 seconds, 730 
nM between 642-702 seconds, and 2200 nM between 760- 
820 seconds. 
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[0106] Before analysing the binding curve, start and end 
times for each injection of analyte are identified, and the 
base line is adjusted to zero level. Generally, the binding 
curve may then be analysed using an interaction model 
based on time dependent analyte concentration and allowing 
for mass transport limitation, such as the above-mentioned 
two-compartment model. If necessary, this model may be 
adapted to other interaction mechanisms, such as, e.g., 
conformation change (A+B<->AB^AB*), competing reac- 
tions (A1+B<-^A1B; A2+B^A2B), parallel reactions 
(A+B1<->AB1; A+B2<-^AB2), or bivalent analyte 
(A+B<-^AB; AB+B<-^AB2). Such models are described in 
more detail in, for example, the BIAevaluation Software 
Handbook mentioned above. 

[0107] Assuming, for example, 1:1 binding, the following 
differential equation may represent the time dependent con- 
centration of analyte A capable of binding to immobilized 
ligand B. 



— = kiiCi + C2 + C3 + ... C„ - A) - {kaA *B- k^AB) 



[0108] where A is the surface concentration of free analyte 
A, is a transport coefficient, to are the different bulk 

concentrations of analyte, B is the concentration of immo- 
bilized ligand (ligand density), AB is the concentration of 
formed complex (=bound analyte), and k^ and k^ are the 
association rate constant and dissociation rate constant, 
respectively. 

[0109] If at the beginning of the experiment the analyte 
concentration. A, at the surface is assumed to be zero, the 
concentration of ligand, B, is assumed to be R^^.^^, and the 
concentration of formed complex, AB, is assumed to be 
zero, the following functions, here written in a form com- 
patible with the above-mentioned BIAevaluation software, 
describe the interaction system. 

AB+RI1*$1+RI2*$2+RI3*$3+RI4*$4+RI5*$5; 

$ 1 ={sign{t-tOnl)-{sign(t-tOffl))/2; 
%2={sign{t-tOn2)-{sign{t-tOffT))l2\ 
$3={sign(t-tOn3)-(sign{t-tOff3))/2; 

$4={sign{t-tOn4)-{sign(t-tOff4y)/2; 
$5=kt'*'{$l'*concl+$2*conc2+$3*conc3+$4*conc4- 

A); 

A=$5-$6\0; 
B=-$6^ax; 

AB=$6\0 

[0110] In the above model, sign functions $1 to $4 (being 
1 or 0) indicate when a certain analyte concentration is 
injected (tOnl is the time when analyte concentration 1 is 
injected, tOffl is the time when analyte concentration 1 is 
stopped, etc). Functions $5 and $6 describe together the 
variation of the surface concentration of analyte over time, 
$6 describing how the ligand density and complex varies 
with time. The function on the first line corresponds to the 
total change in signal from binding AB and bulk effects 
present during each injection. By numerical integration of 
these equations and global fitting to all data points of the 



whole sensorgram in FIG. 3, using, e.g., BIAevaluation 
software (Biacore AB, Uppsala, Sweden) adapted to the 
present invention, the kinetic rate constants k^^ and k^^ can be 
obtained. Mathematical models for the other interaction 
mechanisms may readily be devised by a person skilled in 
the art. 

[0111] The result of the fit is shown in FIG. 4, where the 
solid ''non-rippling line" indicates the adapted curve for best 
fit for the concentration range 80-2200 nM, k^ and k^ being 
constrained to sins le values of l.OnO^ M"^ s"^ and 6.6*10"^ 
s"^, respectively, as described in Example 1 below. 

[0112] With the above model, kinetic and affinity data can 
thus be obtained in spite of transport limitations. Since, as 
mentioned above, the method is insensitive to the injection 
order of analyte, analysis of relatively low analyte concen- 
trations may be performed when binding is close to satura- 
tion or above the equilibrium response, thereby maximizing 
the kinetic information in the sensorgram where transport 
effects are small. 

[0113] Data from binding analyses at two or more different 
ligand densities may advantageously be fitted simulta- 
neously. 

[0114] Measurement of concentration simultaneously with 
kinetic rate constants may be performed, using the above 
titration procedure and data fitting, provided that the trans- 
port coefficient k^ is known and kept constant. Preferably, 
binding data from at least two, and preferably more ligand 
densities are then used. 

[0115] A similar titration type approach as those described 
above may also be used to determine (at least approximate) 
kinetic rate constants for a plurality of different analytes 
binding to the same ligand (at the same or different known 
concentrations). In this case, different analytes are sequen- 
tially passed over the sensor surface rather than different 
analyte concentrations. The evaluation of the binding data is 
modified to permit fitting of different rate constants to the 
resulting total binding curve. Such a procedure may, for 
example, be used to screen different binding partners 
towards the ligand to quickly find the best binder, e.g., drug 
candidates binding to a drug target, such as a receptor. 
Optionally, this approach may be combined with the use of 
different analyte concentrations. 

[0116] As mentioned above, analyte injections may also 
be performed by a "pulse injection" approach. An exemplary 
system for performing such a pulse injection method, based 
on the microfluidic system of a BIACORE® 3000 instru- 
ment (Biacore AB, Uppsala, Sweden), is schematically 
illustrated in FIG. 5. 

[0117] As can be seen in FIG. 5, there are provided two 
vessels (e.g., test tubes) 20 and 21 containing sample (20) 
and buffer (21), respectively. There is also provided a means 
22 for aspirating liquid from the test tubes, indicated with 
vertical lines extending down into the test tubes. This means 
22 can suitably be a needle, and since the same needle is 
used for both liquids, the needle shown in the sample tube 
20 is shown by a broken line. The needle would thereby be 
physically moved between the tubes for the aspiration of 
liquids sequentially. Of course, there are other possibilities 
of devising the aspiration means, the one shown being only 
exemplary. 
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[0118] A system buffer supply is also provided. Initially 
the entire system is filled with buffer, i.e., all tubing contains 
this buffer. Respective segments of tubing 23 (sample and 
system buffer, respectively) are coupled to a microfluidic 
device 24 (in the BIACORE® instruments referred to as an 
Integrated Fluidic Cartridge — ^IFC) enabling controlled liq- 
uid delivery to one or more flow cells, in the illustrated case 
four flow cells 25a to 25d. Each flow cell has a sensor 
surface onto which one or more suitable target(s) are immo- 
bilized. There are also provided a number of valves vl to v4 
in the microfluidic device 24 for the control of the flows of 
the respective liquids. Valves v2 and v3 control the supply 
of fluid to the flow cells whereas valves vl and v4 each are 
connected to waste. The valves vl to v4 are controlled by a 
control unit 26. Alternatively, the flow in the various lines 
can be controlled by accurate pumps, whereby the actual 
flow rates can be monotonically controlled to provide the 
desired flow rates, ranging from zero flow to the maximum 
flow rates required, or combinations thereof. 

[0119] The first step in the procedure is to aspirate a small 
volume of buffer into the needle 22, i.e., to immerse the 
needle into the buffer tube 21, and to aspirate the appropriate 
volume into the needle. It is however not strictly necessary 
to fill the needle with buffer by aspiration. Instead, the 
needle can be filled with buffer from the other end, i.e., from 
the system buffer supply by filling the entire system with 
buffer. Then, the needle is moved to the sample tube and a 
suitable volume of sample of about 500 jul is aspirated. 
However, the actual volume may depend on the application 
and the kind of sample, and can vary within wide limits, say 
between 1 jul and 4 ml. 

[0120] The aspiration of sample will lead to mixing of the 
sample and buffer by dispersion, thereby creating a gradient 
in the tubing. In this case, the gradient will be a decreasing 
gradient (as seen from the needle) running through the 
tubing. If an increasing gradient is required, one would have 
to aspirate buffer after the sample aspiration, and ensure that 
a non-dispersed sample trailing edge is provided by first 
aspirating an air bubble to protect the sample from liquid 
already present in the needle, second a sample and third a 
buffer segment. The aspiration sequence always ends with 
aspiration of one or a few air bubbles to protect the gradient 
from liquid already present in the microfluidic device 24. 

[0121] Prior to the first step, it is preferable to perform a 
few alternating air and sample aspirations to provide con- 
secutive segments of air and sample and to inject them into 

the microfluidic device 24. In this way the sample will be 
protected from unwanted dispersion with running buffer in 
the microfluidic device, i.e., the leading front of the aspi- 
rated sample liquid will exhibit the nominal (maximum) 
concentration. 

[0122] When a gradient has been established, it is injected 
via the needle into the microfluidic device and valves v2 and 
v3 in the sample and buffer flow lines are opened and closed 
according to a programmed sequence to enable alternating 
sample (exhibiting a gradient in the longitudinal direction of 
the tubing) and buffer pulses to be fed into the flow ceUs, 
such that the sample liquid flow is intersected at least once, 
preferably a plurality of times, by a further liquid, repre- 
sented by the system buffer in this case. This intersection 
wiU create at least two separated segments of liquid. How- 



ever, other further liquids than the system buffer are of 
course possible, such as pure solvent, solutions containing 
other species of interest, etc. 

[0123] Thus, the leading edge of a decreasing sample 
gradient flow will represent a first concentration. Most often 

the concentration at the leading edge will be very close to the 
nominal, and can be taken to represent a known concentra- 
tion. However, the major part of the sample flow will exhibit 
a gradient, and thus the majority of said segments that are 
created will have different concentrations with respect to the 
sample. 

[0124] After a predetermined volume of sample gradient 
flow has passed into the flow cells, valve v2 is closed and 
valve v3 is opened, thereby injecting buffer into the line 
behind the sample flow. During the passage of analyte- 
containing sample over the sensor surface having targets 
immobilized on it, the analyte will associate with the targets. 
The volume of sample should preferably be sufficient to 
enable equilibrium to establish. However, it is not always 
required that equilibrium be reached but an equilibrium level 
can be calculated from the corresponding binding curve 
graph (response vs time). The time frames involved depend 
on sample specific binding and transport characteristics, 
flow rate, temperature, flow cell dimensions, etc. 

[0125] When sample has been injected for a sufficiently 

long time, buffer is injected by opening valve v3 and closing 
valve v2. During the passage of buffer over the surface, 
analyte will dissociate. The process is repeated until the 
aspirated sample has been injected. 

[0126] It is not necessary to inject the complete gradient 
into the flow cefl. During buffer injection (v3 open, v2 
closed) valve 1 can be opened to discard a small segment of 
the gradient. This will reduce the number of pulses produced 
and reduce the time needed for a full injection. 

[0127] In an alternative embodiment, wherein a system 
without valves is used, during the association phase, i.e., 
during the time that the sample is passed through the sensor 
cell, the buffer flow is set to a very low value, less than 5%, 
and, e.g., about 1% of the regular flow. This is not strictly 
necessary, but prevents sample solution from leaking into 
the buffer line. Then a certain, predetermined amount of 
sample is injected into the microfluidic device at a specified 
rate. The buffer flow rate is then reset to the regular rate. 
During the passage of buffer through the ceU, sample 
compound that has bound to the target on the sensor surface 
is allowed to dissociate for a suitable time period. 

[0128] In FIG. 6 there is shown an alternative method of 
creating a gradient using the microfludic device in FIG. 5. 
This is done by aspirating a sample segment of known 
concentration and diluting it with buffer in the microfluidic 
device using a connection cl and a tubing segment ml. This 
allows the buffer and sample to form a homogenous mixture 
prior to contacting the flowcells. The pulses would be 
generated as previously described, i.e., by using alternating 
pumps or valves v2 and v3. The connection cl could be a 
simple T-connection so that the concentration of the sample 
in the gradient is controlled by how the ratio of [flowrate 
(buffer)] and [flowrate (sample)] changes over time. Another 
possibility could be to have a two-way valve as connection 
cl. The concentration of the sample will be controlled by 
switching the inlet to ml between buffer and sample, having 
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the two-way valve open for buffer a different time than open 
for sample. In the tubing segment ml the discrete connected 
segments of sample and buffer will form a homogenous mix 
due to dispersion. This method makes it possible to generate 
a gradient with known concentration of the sample at all 
times, in contrast to the dispersion concentration gradient 
where only the first few pulses have a known compound 
concentration. In the latter case, it will be necessary to fit the 
concentrations locally simultaneously as fitting the kinetic 
parameters. This will be described in more detail below. 

[0129] While it currently is preferred to use so-called 
label-free detection methods, including those mentioned 
above, such as SPR, other detection techniques may, as 
already mentioned above, of course, also be contemplated, 
such as those based on the detection of a label, e.g., a 
radiolabel, fluorophore (including surface fluorescence), 
chromophore, chemiluminescer, a marker for scattering 
light, an electrochemically active marker, a thermoactive 
marker, etc. 

[0130] The methods of the invention, including titration 
(i.e., sequential analyte injections) and evaluation with time 
dependent anal3^e concentration, are, as already mentioned 
above, suitably reduced to practice in the form of a computer 
system running software which implements the different 
steps of the procedure. Preferably, the titration is adaptive as 
mentioned above. The invention also extends to computer 
programs, particularly computer programs on or in a carrier, 
adapted for putting the titration procedure of the invention 
into practice. The carrier may be any entity or device capable 
of carrying the program. For example, the carrier may 
comprise a storage medium, such as a ROM, a CD ROM or 
a semiconductor ROM, or a magnetic recording medium, for 
example a floppy disc or a hard disk. The carrier may also 
be a transmissible carrier, such as an electrical or optical 
signal which may be conveyed via electrical or optical cable 
or by radio or other means. Alternatively, the carrier may be 
an integrated circuit in which the program is embedded. 

[0131] Software like that outlined above may also be used 
to generate theoretical interaction curves for different inter- 
action models similarly as described further above for the 
conventional single injection analysis procedures. By gen- 
erating data for a complex model, for instance a conforma- 
tion change model, varying parameters such as injection 
time and analyte concentration and then fitting this data to 
both a 1:1 binding model and a conformational change 
model, the residuals can identify conditions for optimal 
designs of kinetic experiments. 

[0132] In the following Examples, various aspects of the 
present invention are disclosed more specifically for pur- 
poses of iUustration and not limitation. 

EXAMPLES 
[0133] Instrumentation 

[0134] A BIACORE® S51 instrument (Biacore AB, Upp- 
sala, Sweden) was used in Examples 1 to 4 below. This 
instrument has two Y-type flow cells which allow a dual flow 
of fluids over a sensor chip surface, so-called hydrodynamic 
addressing, as described in, for example, EP-B 1-1021703 
(the disclosure of which is incorporated by reference herein). 
The instrument uses three parallel detection spots on the 
sensor chip. 

[0135] A BIACORE®3000 instrument (Biacore AB, Upp- 
sala, Sweden) was used in Example 5 below. In this instru- 



ment, a micro-fluidic system passes samples and running 
buffer through four individually detected flow cells (one by 
one or in series). 

[0136] As sensor chip was used Series S Sensor Chip CMS 
(Biacore AB, Uppsala, Sweden) which has a gold-coated 
surface with a covalently linked carboxymethyl-modified 
dextran polymer hydrogel. 

[0137] The outputs from the instruments via the instru- 
ment control software are "sensorgrams" which are a plots 
of detector response (measured in "resonance units", RU) as 
a function of time. An increase of 1000 RU corresponds to 
an increase of mass on the sensor surface of approximately 
1 ng/mm^. In Examples 1 to 4, evaluation was performed 
using BIAevaluation software, version 3.1 (Biacore AB, 
Uppsala, Sweden), adapted to a kinetic model as described 
above. In Example 5, evaluation was performed using 
BIAevaluation software, version 3.1 (Biacore AB, Uppsala, 
Sweden), Matlab version 5.3 (The MathWorks, Inc. Natick, 
Mass., U.S.A.) and Excel 97 (Microsoft Corp., Redmond, 
Wash., U.S.A.). 

Example 1 

Binding of a Kinase Inhibitor to Immobilized 
Kinase 

[0138] Anti-histidine antibody (Qiagen 34660; Qiagen, 
Venlo, Netherlands) was diluted to 20 jug/ml in 10 mM 
acetate buffer pH 5.0, and 13000 RU of antibody was 
immobilized to a Series S Sensor Chip CM5 using amine 
coupling (Amine Coupling Kit — Biacore AB, Uppsala, 
Sweden) according to the manufacturer's instructions. The 
flow was 10 /il/min. Histidine tagged kinase (proprietary) 
was then captured and cross-linked to the sensor chip 
surface by injecting 30 //g/ml of the kinase in 10 mM PBS 
for three minutes followed by 2.5 min injections of EDC/ 
NHS and ethanolamine to stabilise the kinase on the surface. 
With this procedure 4000 RU of kinase were immobilized. 

[0139] A kinase inhibitor, staurosporin, was injected at 
concentrations of 80, 240, 730 and 2200 nM at a flow of 30 
/il/min. The assay buffer was 50 mM Tris pH 7.5, 150 mM 
NaCl, 10 mM MgCl2 and 3% DMSO. Overlay plots of 
observed and fitted data for the injection sequence are 
demonstrated in FIG. 3 (observed) and FIG. 4 (observed 
plus fitted). The association rate constant was calculated to 
1,0*10^ M"^ s~^, and the dissociation rate constant k^ to 
6,6*10 s-\ 

Example 2 

Binding of a Set of Carbonic Anhydrase Inhibitors 
to Immobilized Carbonic Anhydrase 

[0140] Carbonic anhydrase (Sigma C3640), 20 jug/m\ in 
10 mM acetate buffer pH 5.0, was immobilized to a Series 
S Sensor Chip CM5 using amine coupling according to the 
manufacturer's instructions (flow 10 //1/min), to obtain sen- 
sor surfaces with approximately 2000 and 700 RU, respec- 
tively, of immobilized carbonic anhydrase. Kinetic analysis 
was performed by injecting inhibitors in 10 mM PBS buffer 
with 0.005% P20 and 5% DMSO at a flow of 30 /il/min. The 
inhibitors were acetazoleamide (Sigma A6011), azosulfa- 
mide (Sigma A2759), and benzenesulfonamide (Aldrich 
10,814-6). Overlay plots demonstrating observed and fitted 
data for the two ligand densities are presented in FIG. 7 
(acetazoleamide), FIG. 8 (azosulf amide) and FIG. 9 (ben- 
zenesuffonamide). In FIG. 7, data represent two repeats on 
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each ligand density, and was calculated to 1,8*10^ M ^ 
s-\ and kd to 4,0*10"^ s-\ In FIG. 8, was calculated to 
2,9*10^ s-\ and k^ to 8,1*10"^ s-\ and in FIG. 9 k, 
was calculated to 1,6*10^ M"^ s-\ and k^ to 1,4*10"^ s-\ 

Example 3 

Binding of a High Affinity Inhibitor to Immobilized 
Thrombin 

[0141] Thrombin (Sigma T 1063), 20 /<g/ml in 10 mM 
acetate buffer pH 5.0, was immobilized to a Series S Sensor 
Chip CMS using amine coupling according to the manufac- 
turer's instructions (flow 10 /il/min), resulting in sensor 
surfaces with approximately 2300 and 600 RU, respectively, 
of immobilized thrombin. Kinetic analysis was performed 
by injecting a high affinity thrombin inhibitor, melagatran (a 
gift from Johanna Deinum, AstraZeneca, Molndal, Sweden), 
in 10 mM PBS buffer with 0.005% P20, 5% DMSO and 3.4 
mM EDTA at a flow of 30 /A/min. Overlay plots demon- 
strating observed and fitted data are presented in FIG. 10. 
The calculated k^ was 2,0n0^ M"' s-\ and k^ 1,8*10"^ s-\ 

[0142] Examples 1 to 3 above demonstrate that a wide 
range of rate constants can be determined with the sequential 
injection procedure. In these Examples k.^ values range from 
2,9* 10^ to 2,0* 10^ M"^ s"^ and k^ values range from 0,14 to 
6,6 lO"'^ s""^. The Examples also illustrate that kinetic data 
can be obtained using one or two ligand densities and that 
the injection order can be mixed and not only in the order 
from low to high concentration. 

Example 4 

Combined Sequential and Single Injections 

[0143] In FIG. 11, simulated data representing sequential 
and single injections are combined in the evaluation where 
all binding curves are fitted simultaneously. The figure 
shows an overlay plot of four binding curves and fitted data. 
The three curves to the left represent sensorgrams with 
single injections of analyte at 51200 nM but with varying 



injection times, 10 s, 30 s and 500 s, whereas the curve to 
the right represents sequential injections of analyte from 800 
to 12800 nM. This illustrates that sensorgrams with varying 
numbers of injections and injection times can be analysed. 

Example 5 

Estimation of Interaction Rate Constants 

[0144] The kinetics of a camel derived heavy chain triple 
mutant single domain antibody (cAb-Lys3:s SGS) binding 
to lysozyme (obtained from the Department of Ultrastruc- 
ture, Vrije Universiteit, Brussels, Belgium) was studied with 
the "pulse injection method" described above. All experi- 
ments were performed at 30° C. 190 RU (chip 1) and 280 
RU (chip 2) of lysozyme was immobilized using a standard 
amine coupling procedure. Upon a 2-min activation, 
lysozyme (8 jug/ml in 10 mM Na2HP04, pH 7.0) was 
injected for 3 min (chip 1) and 4 min 30 s (chip 2). 10 mM 
Na2HP04, pH 7.0 (flow rate 5 juVmin) was used as running 
buffer during immobilization. SGS was injected at different 
initial concentrations (0.5, 1.0 and 2.0 juM in HBS-EP). The 
results are shown in FIG. 12. 

[0145] Bulk errors in the sample solutions were corrected 
for by subtraction of the reference flow cefl signals. Indi- 
vidual pulses were separated and aligned, using MATLAB, 
so that each pulse corresponded to one binding curve. The 
curves were superimposed in the BIAevaluation software. 
15 pulses were used in every fit. The first two pulses were 
assumed to be of initial concentration. Global starting values 
of k^, kd and R^ax were fitted to the second pulse (pulse 
number one was omitted because of its irregular shape), 
since its concentration was known. These values were then 
used to locally fit the concentrations of all pulses, k^, k^ and 
Rmax estimations were refined, using the new concentration 
information. The process was repeated until all parameters 
converged. Each pulse injection was evaluated separately. 
The fitting of the concentration resulted in a partially linear 
concentration gradient. Kinetic data obtained with the pulse 
injection method is presented together with mean values and 
standard deviations in Table 1 below. 
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*Co is the nominal concentration of SGS 

* is a statistical measure of the quality of the fit 
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[0146] It is to be understood that the invention is not 
limited to the particular embodiments of the invention 

described above, but the scope of the invention will be 
established by the appended claims. 

1. A method of determining kinetic parameters for a 
reversible molecular interaction between a ligand immobi- 
lized to a solid support surface and a binding partner to the 
ligand in solution, comprising the steps of: 

a) sequentially, without intermediate regeneration or 
renewal of the immobilized ligand, flowing a plurality 
of fluid volumes containing different known concen- 
trations of the binding partner over the solid support 
surface to permit association of binding partner to the 
immobilized ligand, 

b) flowing over the solid support surface a fluid volume 
free from binding partner to permit dissociation of 
binding partner from the ligand, 

c) monitoring during steps a) and b) the momentary 
amount of binding partner bound to the solid support 
surface related to time and solution concentration of 
binding partner and collecting the binding data, and 

d) determining the kinetic parameters by fitting a prede- 
termined kinetic model for the interaction between the 
binding partner and the immobilized ligand to the 
collected binding data, which model allows for mass 
transport limitation at the solid support surface. 

2. The method according to claim 1, wherein step d) 
comprises determining the kinetic parameters by globally 
fitting the kinetic model to the collected binding data. 

3. The method according to claim 1, wherein steps a) to 
c) are repeated at least once with the same or another 
sequence of binding partner concentrations, and all the 
collected binding data are fitted simultaneously to the kinetic 
model in step d). 

4. The method according to claim 1, wherein sequential 
binding data obtained from steps a) to c) are combined with 
at least one set of binding data for a single concentration of 
binding partner, and the combined binding data are fitted in 
step d) to the kinetic model. 

5. The method according to claim 1, wherein step b) is 
performed after flowing the last fluid volume containing 
binding partner over the surface. 

6. The method according to claim 1, wherein step b) is 
performed between two fluid volumes containing binding 
partner. 

7. The method according to claim 1, wherein each volume 
containing binding partner is followed by a fluid volume free 
from binding partner. 

8. The method according to claim 7, wherein one of the 
fluid volumes free from binding partner is flowed over the 
solid support surface for a substantially longer time than the 
other fluid volumes free from binding partner. 

9. The method according to claim 8, wherein the last fluid 
volume free from binding partner is flowed over the solid 
support surface for a substantially longer time than the other 
fluid volumes free from binding partner. 

10. The method according to claim 1, wherein the plu- 
rality of fluid volumes containing binding partner include 
fluid volumes with the same concentration of binding part- 
ner. 



11. The method according to claim 1, wherein the fluid 
volumes containing binding partner contain successively 
increasing concentrations of binding partner. 

12. The method according to claim 1, wherein step a) is 
stopped when an increased concentration of binding partner 
results in a predetermined marginal increase in the amount 
of bound binding partner at the sensor surface. 

13. The method according to claim 1, wherein the kinetic 
model is a two-compartment model. 

14. The method according to claim 1, wherein the kinetic 
model is adapted to an interaction mechanism selected from 
1:1 binding, parallel reactions, competing reactions, confor- 
mation change and bivalent analyte. 

15. The method according to claim 1, wherein step a) is 
performed at at least one other density of ligand at the solid 
support surface. 

16. The method according to claim 1, wherein the flow 
over the solid support surface is maintained essentially 
constant. 

17. The method according to claim 1, wherein the fluid 
volumes containing different known concentrations of the 
binding partner comprise volume segments taken from a 
concentration gradient of the binding partner. 

18. The method according to claim 1, wherein the solid 
support surface is a sensor surface. 

19. The method according to claim 1, wherein the detec- 
tion at the solid support surface is based on evanescent wave 
sensing. 

20. The method according to claim 19, wherein the 
detection at the solid support surface is based on surface 
plasmon resonance (SPR). 

21. The method according to claim 1, wherein the detec- 
tion at the solid support surface is based on detection of a 
label. 

22. The method according to claim 1, wherein step d) is 
computer-implemented. 

23. The method according to claim 1, wherein steps a) to 
d) are computer- implemented. 

24. A method of determining the concentrations of bind- 
ing agent in a plurality of samples, wherein the binding 
agent is capable of interacting with a ligand immobihzed to 
a solid support surface, comprising the steps of: 

a) sequentially, without intermediate regeneration or 
renewal of the immobilized Hgand, flowing a plurality 
of different sample volumes containing binding agent 
over the solid support surface, 

b) flowing over the solid support surface having binding 
agent bound thereto a fluid volume free from binding 
agent to permit dissociation of binding agent from the 
ligand, 

c) monitoring during steps a) and b) the momentary 
amount of binding agent bound to the solid support 
surface related to time and solution concentration of 
binding agent and coUecting the binding data, and 

d) determining the concentrations of binding agent in the 
different sample volumes by fitting a predetermined 
kinetic model for the interaction between the binding 
agent and the immobilized ligand to the collected 
binding data, which model allows for mass transport 
limitation at the solid support surface, and wherein the 
mass transport coefficient is known and kept constant. 
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25. The method according to claim 24, wherein step d) 
comprises determining the concentrations by globally fitting 
the kinetic model to the collected binding data. 

26. The method according to claim 24, wherein also 
kinetic parameters are determined. 

27. The method according to claim 24, wherein the solid 
support surface is a sensor surface. 

28. The method according to claim 24, wherein the 
detection at the solid support surface is based on evanescent 
wave sensing. 

29. The method according to claim 28, wherein the 
detection at the solid support surface is based on surface 
plasmon resonance (SPR). 

30. The method according to claim 24, wherein the 
detection at the solid support surface is based on detection 
of a label. 

31. The method according to claim 24, wherein step d) is 
computer-implemented. 

32. The method according to claim 24, wherein steps a) to 
d) are computer-implemented. 

33. A method of determining kinetic parameters for 
reversible molecular interactions between a ligand immobi- 
lized to a solid support surface and a plurality of different 
binding partners to the ligand in solution, comprising the 
steps of: 

a) sequentially, without intermediate regeneration or 
renewal of the immobilized ligand, flowing a plurality 
of fluid volumes containing different binding partners 

over the solid support surface to permit association of 
binding partners to the immobilized ligand, 

b) after the addition of each fluid volume containing 
binding partner, flowing over the solid support surface 
a fluid volume free from binding partner to permit 
dissociation of binding partner from the ligand, 

c) monitoring during steps a) and b) the momentary 
amount of binding partner bound to the solid support 
surface related to time and solution concentration of 
binding partner and collecting the binding data, and 

d) determining the kinetic parameters by fitting a prede- 
termined kinetic model for the interaction between each 
binding partner and the immobilized ligand to the 
collected binding data, which model allows for mass 
transport limitation at the solid support surface. 

34. The method according to claim 33, wherein step d) 
comprises determining the kinetic parameters by globally 
fitting the kinetic model to the collected binding data. 

35. The method according to claim 33, wherein the solid 
support surface is a sensor surface. 

36. The method according to claim 33, wherein the 
detection at the solid support surface is based on evanescent 
wave sensing. 

37. The method according to claim 36, wherein the 
detection at the solid support surface is based on surface 
plasmon resonance (SPR). 

38. The method according to claim 33, wherein the 
detection at the solid support surface is based on detection 
of a label. 

39. The method according to claim 33, wherein step d) is 
computer-implemented. 

40. The method according to claim 33, wherein steps a) to 
d) are computer-implemented. 



41. A method of determining kinetic parameters for a 
reversible molecular interaction between a Hgand immobi- 
lized to a sensor surface and a binding partner to the ligand 
in solution, comprising the steps of: 

a) sequentially, without intermediate regeneration or 
renewal of the immobilized ligand, flowing a plurality 
of fluid volumes containing different concentrations of 
the binding partner over the solid support surface to 
permit association of binding partner to the immobi- 
lized ligand, wherein the fluid volumes containing 
binding partner are discrete segments of a concentra- 
tion gradient of the binding partner separated by seg- 
ments of fluid free from binding partner, 

b) monitoring during step a) the momentary amoimt of 
binding partner bound to the solid support surface 
related to time and solution concentration of binding 
partner and collecting the binding data, and 

c) determining the kinetic parameters by fitting a prede- 
termined kinetic model for the interaction between the 
binding partner and the immobilized ligand to the 

collected binding data. 

42. The method according to claim 41, wherein the kinetic 
model allows for mass transport limitation at the solid 
support surface. 

43. The method according to claim 41, wherein the 
binding partner concentrations of all gradient segments are 
known, and a global fitting is performed in step c). 

44. The method according to claim 41, wherein the 
binding partner concentrations of all gradient segments are 
not known, and the fitting performed in step c) includes local 
fitting of the binding partner concentration. 

45. The method according to claim 41, wherein the solid 
support surface is a sensor surface. 

46. The method according to claim 41, wherein the 
detection at the solid support surface is based on evanescent 
wave sensing. 

47. The method according to claim 46, wherein the 

detection at the solid support surface is based on surface 
plasmon resonance (SPR). 

48. The method according to claim 41, wherein the 
detection at the solid support surface is based on detection 
of a label. 

49. The method according to claim 41, wherein step c) is 
computer-implemented. 

50. The method according to claim 41, wherein steps a) to 
c) are computer-implemented. 

51. An analytical system for detecting molecular binding 
interactions, comprising: 

(i) a sensor device comprising at least one sensor surface, 
detection means for detecting molecular interactions at 
the at least one sensor surface, and means for producing 
detection data representing the progress of the interac- 
tions with time, and 

(ii) data processing means for performing the steps as 
defined in claim 1. 

52. An analytical system for detecting molecular binding 

interactions, comprising: 

(i) a sensor device comprising at least one sensor surface, 
detection means for detecting molecular interactions at 
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the at least one sensor surface, and means for producing 
detection data representing the progress of the interac- 
tions with time, and 

(ii) data processing means for performing the steps as 
defined in claim 24. 

53. An analytical system for detecting molecular binding 
interactions, comprising: 

(i) a sensor device comprising at least one sensor surface, 
detection means for detecting molecular interactions at 
the at least one sensor surface, and means for producing 
detection data representing the progress of the interac- 
tions with time, and 

(ii) data processing means for performing the steps as 
defined in claim 33. 

54. An analytical system for detecting molecular binding 
interactions, comprising: 

(i) a sensor device comprising at least one sensor surface, 
detection means for detecting molecular interactions at 

the at least one sensor surface, and means for producing 
detection data representing the progress of the interac- 
tions with time, and 

(ii) data processing means for performing the steps as 
defined in claim 41. 

55. A computer program comprising program code means 
for performing the method of claims 1 when the program is 
run on a computer. 

56. A computer program comprising program code means 
for performing the method of claims 24 when the program 
is run on a computer. 

57. A computer program comprising program code means 
for performing the method of claims 33 when the program 
is run on a computer. 

58. A computer program comprising program code means 
for performing the method of claims 41 when the program 
is run on a computer. 



59. A computer program product comprising program 
code means stored on a computer readable medium or 
carried on an electrical or optical signal for performing the 
method of claim 1 when the program is run on a computer. 

60. A computer program product comprising program 
code means stored on a computer readable medium or 
carried on an electrical or optical signal for performing the 
method of claim 24 when the program is run on a computer. 

61. A computer program product comprising program 
code means stored on a computer readable medium or 
carried on an electrical or optical signal for performing the 
method of claim 33 when the program is run on a computer. 

62. A computer program product comprising program 
code means stored on a computer readable medium or 
carried on an electrical or optical signal for performing the 
method of claim 41 when the program is run on a computer. 

63. A computer system containing a program for perform- 
ing the method of any one of claim 1. 

64. A computer system containing a program for perform- 
ing the method of any one of claim 24. 

65. A computer system containing a program for perform- 
ing the method of any one of claim 33. 

66. A computer system containing a program for perform- 
ing the method of claim 41. 

67. A method for simulating theoretical molecular inter- 
action data, wherein a computer program according to claim 

55 is used. 

68. A method for simulating theoretical molecular inter- 
action data, wherein a computer program according to claim 

56 is used. 

69. A method for simulating theoretical molecular inter- 
action data, wherein a computer program according to claim 

57 is used. 

70. A method for simulating theoretical molecular inter- 
action data, wherein a computer program according to claim 

58 is used. 



